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Abstract 

Topological quantum phase transitions intrinsically intertwine self-similarity and topology of 
many-electron wave-functions, and divining them is one of the most significant ways to advance 
understanding in condensed matter physics. Our focus is to investigate an unconventional class of 
the transitions between insulators and Dirac semimetals whose description is beyond conventional 
pseudo relativistic Dirac Hamiltonian. At the transition without the long-range Coulomb interac¬ 
tion, the electronic energy dispersion along one direction behaves like a relativistic particle, linear 
in momentum, but along the other direction it behaves like a non-relativistic particle, quadratic in 
momentum. Various physical systems ranging from Ti02-V02 heterostructure to organic material 
a— (BEDT-TTF) 2 l 3 under pressure have been proposed to have such anisotropic dispersion rela¬ 
tion. Here, we discover a novel quantum criticality at the phase transition by incorporating the 
1 long range Coulomb interaction. Unique interplay between the Coulomb interaction and elec¬ 
tronic critical modes enforces not only the anisotropic renormalization of the Coulomb interaction 
but also marginally modified electronic excitation. In connection with experiments, we investigate 
several striking effects in physical observables of our novel criticality. 
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Quantum criticality and topology are two of the main impetuses of modern condensed 
matter physics. Self-similarity of many-electron wave-functions associated with quantum 
criticality unveils emergent universality of physical observables, and topology of the 
electronic wave-functions manifests itself as various fascinating topological insulators and 
associated quantized responses [1H7]. The two striking characteristics of the wave-function 
are naturally and inevitably intertwined at topological quantum phase transitions. 

Long-range 2 Coulomb interaction between electrons induces striking screening effects 
near the topological phase transitions. Electronic critical modes and the Coulomb interaction 
are intrinsically correlated , so non-trivial quantum criticality usually appears [8ffl3]. For 
example, quasi-particles lose their stability due to the Coulomb interaction and the ground 
state becomes quantum critical non-Fermi liquid with emergent full rotational symmetry 
in quadratic band touching semimetals, which is near three dimensional (3d) topological 
insulator [8]. 

In two dimensions (2d), the Coulomb interaction becomes more special. It is because the 
Coulomb potential originally lives in 3d but electrons are confined in 2d. Thus the electrons 
in 2d feels the dimensionally different interaction, originating from 3d. Since correlation and 
fluctuation are enhanced in lower dimensions, one may expect stronger interplay between 
the Coulomb interaction and critical modes in a topological phase transition, and indeed we 
find the novel quantum criticality in a class of 2d topological quantum phase transitions. 

Conventional 2d topological phase transitions between two topologically distinct insu¬ 
lators are described by the pseudo-relativistic Dirac fermion theory H D (k) = v x k x a x + 
v y k y a y + Ma z , with Pauli matrices in band index spinor space. Here the topological na¬ 
ture of the transition is captured by the change in the Berry curvature of the wave-function 
depending on the sign of M, and different patterns of opening up band gaps at separate 
Dirac points represent different topological insulator phases when supplemented with proper 
symmetries. The long-range Coulomb interaction at the critical point (M = 0) induces in¬ 
triguing logarithmic modification of the Dirac velocities, so not only rotational symmetry at 
the critical point emerges but also important interaction effects appear whose structure has 
been extensively studied in literature mm in connection with charge-neutral mono-layer 
graphenes. We emphasize that the isotropic 2 Coulomb interaction dominates microscopic 
anisotropy of electrons in this case. 

Here we focus on a different class of the topological phase transitions whose electronic 
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Hamiltonian is 


H(k) = vk x a x + Ak 2 y a y . (1) 

With a tuning parameter m, energy spectrum of H = H(k)—mo y is E±(k ) = ±^jv 2 k 2 + (Ak 2 — m ) 2 . 
The two phases are determined by the sign of the tuning parameter. With a positve A and 
m < 0, energy spectrum is gapped, so the ground state is an insulator. On the other 
hand, with m > 0 , the zero-energy points appear at two points in momentum space, 
k y = ±k* = ±/f. By expanding the Hamiltonian ([TJ) near these points, we obtain 
H + (p) = vp x a x + Vyp y cr y around the point k = ( p x ,k* + p y ) with v y = 2 Ak*. Thus the 
ground state is a 2d Dirac semimetal. Thus, it is clear that our model Hamiltonian describes 
a phase transition between a (either topological or trivial) insulator and a Dirac semimetal 
in 2 d. 

The Hamiltonian ([TJ) has been suggested in various physical systems, ranging from Ti 02 - 
VO 2 oxide heterostructures [TEHTS] and the organic material a — (BEDT-TTF) 2 1 3 under 
pressure [I9H2T] to optical lattice systems j 22 i [24] . For example, in the oxide heterostructure 
Ti 02 -V 02 layers [T 6 HT 8 ] . there is a metal-insulator transition as the number of layers is 
changed. At the certain number of layers, the first-principle band structure calculation[16F 
IT 8 ] reveals that there should be the anisotropic semimetal ([Tj). Furthermore, the structure 
of the Hamiltonian ([TJ) is similar to that of the notorious quantum criticality problem with 
Fermi surfaces in 2 d [25H3D] whose scaling of the dispersion along the radial direction to 
the Fermi surface is linear in momentum while that along the perpendicular direction is 
quadratic in momentum. Due to the similarity in dispersions, we expect that our analysis 
might shed some light on understanding the quantum criticality with Fermi surfaces despite 
of the finite density of states in the Fermi surface case. 

In this work, we show, by using the systematic renormalization group (RG) method, 
that the long range interaction strongly changes the nature of the eigenstates of the 
non-interacting Hamiltonian. We find a novel quantum criticality characterized by both 
anisotropically renormalized and marginal Coulomb interactions which is in sharp contrast 
to other quantum criticalities. The anisotropic marginal quantum criticality is out of intri¬ 
cate interplay between the long range Coulomb interaction in 2d and the critical electron 
modes, and we emphasize its striking properties by calculating physical observables. 
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MODELS WITH COULOMB INTERACTION 


We start with the theory incorporating the electron Hamiltonian with the long-range } 
Coulomb interaction, 


<S — <S^p + <S(/, ( 2 ) 

= j d 2 xdmp\{d r + iecp) + H(—iS7))ip 

S<t> = J d 2 xdr ^0 vT^ 0 = I ^I<7| !</>(<?, ^) I 2 , 

where 0 mediates Coulomb interaction between electrons ip. The short-hand writing J = 
f C (2wP use d- Hereafter, all integrations are defined with the short-distance (or high- 
energy) ultra-violet (UV) cutoff. H(—iV) is the Fourier transformation of H(k) ([Tj), and the 
bare gauge boson propagator g^l{q) — \q\ represents the long-range - Coulomb interaction. 

For future convenience, we introduce a dimensionless coupling constant, the fine structure 
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constant a = which measures the “strength” of Coulomb interaction. 

We investigate the stability of the theory by the lowest order perturbation calculation, 
in particular, by calculating the bosonic self-energy whose Feynman diagram representation 
is in Fig. [I] (a), 


n (g,n) = 


Tr 


' k.LO 


g f {k + q,uj + £l)g f (k, u>) 


gj 1 {k,uj) = —iu> + H{k ) is used. It is straightforward to evaluate the integral (see the 
supplementary information I for detail), and we find that 


!%) - 2 \%\ G(£ q ), P Q - \ ( 3 ) 

where G(£ q ) is the function of the dimensionless parameter t; q . Hereafter, we drop the fre¬ 
quency dependence in the boson self-energy since we are only interested in the instantaneous 
Coulomb interaction. The full functional form of G(£) is not important. Thus we will not 
present it here and plot it only in the supplementary information I. Instead, the asymptotic 
behavior of n(qr) in each direction is extracted 

x olc x v\q x \ . ac v 

n (g*,0) = — —y —, n (0.9w) = — 2 ~ x \q y \ 


with the numeric constants {c x ~ 2.7, c y ~ 2.5}. Notice that the boson self-energy is 
independent of the UV cutoff which signal a novel quantum criticality in our system. 
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FIG. 1 . Diagrams for (a) the boson self-energy n(g) and (b) the fermion self-energy 

Here the dotted line represents the boson propagator and the solid line represents the 

fermion propagator gf(k,uj). 

It is clear that the perturbation becomes more important than the original bare term 
along q = q x x , 


kb,o(9*>°)l ~ \Qx\ < |n(g*,0)| ~ vB) 

in the limit q —>■ 0. Thus we conclude that the action (j2j) is unstable under the fermion-gauge 
boson coupling. 

The instability from the perturbative calculation often indicates the presence of the stable 
strong-coupling fixed point which can be accessed by large- Nf analysis with the number of 
fermion flavors Nf. The large- Nf analysis starts with adding the bosonic self-energy to the 
boson bare term, 

S* ~ NfU^q^q,^ 2 . 

J q,o; ^ 

The schematic representation of the inverse of the corrected boson propagator g^ l {q) = 
\q\ — NfU(q) is 

~ \<l\ + ^{ c y\%\+ c x\J^j^)- ( 4 ) 

The limit a —> 0 recovers the unstable bare action (J2]), and we investigate the opposite limit 
Nfa —> oo where we drop the bare term (~ |q|). 

Using this corrected boson propagator, we calculate the fermion self-energy in Fig. 0(0 

— oj) = Sv( A, fi)k x a x + 5H(A, fi)k 2 y a y , 
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obtained by expanding the self-energy near k = 0 with the UV and IR cutoffs, A, /j. Straight¬ 
forward calculation gives 


v f p A JMf p 

with the two dimensionless constants, J x ~ 0.18, J y ~ 0.03 (see supplementary information 
III for detail). We notice that the instantaneous nature of the Coulomb propagator enforces 
no vertex correction through the Ward identity. 

Therefore, the RG flow equations, i.e., beta functions, for v and a can be derived by 
changing the ratio, - = e l , 


dv 2 J x da 2 J x 

~dl = l^ v ' = 


(5) 


near the strong-coupling fixed point. It is clear from the RG equations that the fine structure 
constant a decreases with the anomalous dimension of the velocity, This concludes that 
the strong-coupling fixed point is unstable. 

Our controlled analysis near the two extreme limits (standard perturbation and large -Nf 
analysis) clearly shows that both the fixed points are unstable. Then it is obvious that the 
stable fixed point should be in the intermediate regime, which is difficult to access in a fully 
controlled way. Thus, we study the fixed point with the standard momentum-shell RG and 
check a posteriori its validity by self-consistency. 

In the momentum-shell RG analysis, we remark that the non-analytic dependence |q| of 
the Coulomb interaction does not receive correction from integrating out higher-momentum 
modes. Thus, we first keep the seemingly irrelevant ~ kq 2 term in the boson action, 




'q,u 


-{\q\ + kqlmq^)\\ 


It turns out that the following three dimensionless parameters determine the RG flows 

e 2 kA 2 A 3 


a = 


7 = 


0 = f ■■ 

37 


( 6 ) 


2tt 2 v ’ v 

Evaluating Feynman diagrams in Fig. [l] gives the renormalized action S' = S' tn + S^. 

Here, we use the cutoff scheme such that we integrate along the ^-directional momentum 

(Ae^ z ,A) with l 1 after integrating out the ^-direction momentum and frequency. On 

integrating out the higher-momentum modes, three parameters are renormalized as 

8v 5 A 5k 

— = aFi 7 l, — = aF 2 7 l, — = fil 
V A K 
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We find that the functions Fi^[ 7 ] are non-negative near 7 = 0 whose specific forms are 
illustrated in supplementary information IV. The RG flow equations of {a, 7 ,/ 3 } are 

da ~ r , 

n = 

‘jj- = 7(0 - 3 + 20(^217] - UW)). 

^=/3(3-^ + a(F 1 [7]-2F 2 [7])). (7) 


The two fixed points are, (a, 7 ,/?) = (0, 0, 0), (0, 0, 3), and it is easy to show that the 
former is unstable and the latter is stable. At the stable fixed point (a, 7 , /3) = (0,0,3), the 
boson propagator receives a large anomalous scaling dimension, which can be understood 
as (1 + f3l) qi ~ q 2 e 31 —> W\q x \ (because e~ 21 is the scaling factor of q x ). Such large 

/ 3=3 

anomalous dimension indicates that the momentum-shell RG is not controlled and a priori 
not reliable. 

At the stable fixed point, the effective bosonic action becomes 


S l = l k + K VVh\)\(l>(<Lu)\ 2 ( 8 ) 

which is very similar to the large -Nj calculation with one important difference; a new cou¬ 
pling constant (k) with UV cutoff scale naturally enters in contrast to the larg e-Nf calcula¬ 
tion where the coefficient of \j'\q x \ is ~ (see equation Q) which depends on the other 
parameters {v, A, a}. Here the new dimensionful parameter n appears in the bosonic part 
at the intermediate coupling regime. 

With this intuition in hand, we investigate the stability of the new fixed point by taking 
equation ([ 8 ]) as the bare boson action and performing the momentum-shell RG near this 
fixed point. Remarkably, we find that the velocity v and inverse mass A receives the same 
corrections at the fixed point 

Sv A 5 A A 

— = a-C f log(-), -=a-C f log(-), 

U il fJj 

The same correction is another evidence for our fixed point to be stable since the ratio 
appearing in the boson self-energy IT(qr) becomes constant. Notice that the remarkable same 
correction also appears in Fermi surface quantum criticalities with very different physical 
reasons which also supports stability of our fixed point [26]. It is manifest that the gapless 
excitation structure of our system is completely different from that of Fermi surfaces (lines) 
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in 2d. Only nodal point excitation appears in our system. However, the low energy scaling 
structures of the two systems are same considering the patch theory of Fermi surfaces in 2d 
; one momentum direction has linear scaling while the other one has quadratic scaling. We 
believe this unexpected similarity is the source of the similar behaviors in the beta functions. 
It would be very intriguing to find more similarity and difference of two systems’ quantum 
criticalities, which we leave for future work. 

The beta functions around the novel fixed point are 


da 

~di 


- C f a 2 , 


dv 

~dl 


a ■ CfV , 


dA 

Id 


a ■ CfA 


(9) 


where Cf ~ 0.8 calculated in the supplementary information V. 

We remark that the hard momentum cutoff scheme is only used for simplicity and illus¬ 
tration. It is shown that our results are independent of cutoff schemes in the supplementary 
information VIII. 

The above RG flow structure ([9]) is unique to this fixed point. The beta functions contain 
the fine structure constant a in contrast to those of the large -Nf calculation ([5]) in which a is 
absent, and here both v, A receive the same logarithmic corrections, which are proportional 
to a. Thus, the fine structure constant a decreases and the fermion only receives the 
logarithmic corrections, which indicates the fixed point is stable. Naturally, as in mono-layer 
graphene, marginal Fermi liquid behaviors are expected with higher order corrections. ra 

Based on the calculations and intuitions, A schematic RG flow can be deduced as in Fig. 
[2] which summarizes our main results. Our controlled calculation shows the non-interacting 
critical point (Non-Int.) and the strong coupled fixed point (S) are unstable and the RG 
flow comes out of the both points and flow into the intermediate fixed point (QC), which 
is characterized by the definite anisotropic scaling of bosons and electrons and the single 
logarithmic corrections to velocity and inverse mass. 


EXPERIMENTAL SIGNATURES 

We now investigate the physical consequences of both the anisotropy and marginal irrel¬ 
evance of the renormalized Coulomb potential at the novel intermediate critical point. 

First of all, with the beta functions of v and A ©. we can find logarithmic corrections 



Ins. non-Int. Dirac 


FIG. 2. Proposed RG flow. The horizontal axis is for the tuning parameter m of the quantum 
criticality equation 0 and the vertical axis is for the strength of Coulomb interaction. There are 
two stable fixed points, insulators (‘Ins’) and Dirac semimetal (‘Dirac’). The two unstable critical 
points are illustrated with dashed circles, non-interacting (‘Non-Int.’) and strong-coupling fixed 
point (‘S’). And the stable critical point is the filled circle (‘QC’). The critical point is characterized 
by the definite anisotropic scaling and the logarithmic corrections to mass and velocity. Near the 
fixed points (Ins, Dirac, QC), one-particle spectrum with the Coulomb interaction is illustrated. 


to all physical quantities. The parameters [v, A} at the temperature scale T are 


v<X) =fo(l + oC/log(h!)). MT) — A>(1 + otC'i log(hi)). 


( 10 ) 


Here E 0 is the bandwidth or the UV cutoff of the theory 0. v 0 and A 0 are the bare 
parameters at the highest energy scale ~ Eq. The logarithmic corrections in {v,A} may be 
observed in quasi-particle experiments as in graphene, for example, angle resolved photo¬ 
emission spectroscopy (ARPES)[32j or quantum oscillation. 

Furthermore, thermodynamic quantities such as specific heat and compressibility also 
show logarithmic corrections. Specific heat and compressibility of the unstable free electron 
fixed point are C V (T) = and = ~ in which E ( T ) 

is the (thermal-averaged) energy density per volume as the function of temperature T and 
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n(/i, T) is the density of the electron per volume as the function of chemical potential /a and 
temperature T. But at the novel fixed point, the logarithmic corrections give 


C V (T) 


T 3 / 2 0.38 

voVAo (1 + otC f log(f?)) 3 / 2 ’ 


k(T) 


0 . 07 ^ 

vqVA){1 + aCf log(f?)) 3 / 2 ’ 


( 11 ) 


at the temperature T by following the reference |[15J. 

Secondly, we can see the effect of the anisotropic renormalization of the gauge boson via 
the screening charge when a single impurity charge Z is introduced at r — 0 . At the level 
of the linear response theory, the screening charge is Pmd(<z) = ZD(q)U 0 (q), in which D(q) 
is the propagator of the gauge boson. We are interested in the directional behaviors of the 
screening charge and hence define the integrated screening charges Q x (x ) = dyp inc i(r) 
and Q y (y ) = f^° dx p in d(n) along x and y. Here we will contrast the extremely different 
behaviors of the screening charges between the free fixed point and the non-trivial fixed 
point. 


At the free fixed point, we ignore the corrections to the gauge boson propagator and use 

Dq( q) = A. Following the straightforward calculation in supplementary information VI, we 

find Q x (x ) oc — j=, Q y (y) oc —Z5(y) whose sign is the opposite of the impurity charge Z. 
v M 


On the other hand, at the non-trivial anisotropic fixed point in which we use the renor¬ 
malized boson propagator D rem {q)~ l — |qr| + K\J\q x \, we find the asymptotic behaviors of 
the screening charges 


Q x (x) oc 


Z 

x|(l + aCf log(|a:/ro|)) 2 ’ 


IsKl + aCaogd^/rol)) 2 ' 


( 12 ) 


where r 0 1 ~ E 0 is the UV cutoff. Here the sign of the screening charge is the same as the 
impurity charge Z , which is reminiscent of graphene case|33j. 

From the above calculations, we see that the asymptotic scaling behaviors of the screening 
charges in distance from the impurity along x and y are surprisingly isotropic. The isotropic 
scaling in both the directions is originated from the facts that the scaling of D rem (q) is 
identical to that of n 0 (q) and that {n(/i), A(/r)} at the energy scale receive the same 
logarithmic corrections as (§• Hence this isotropic scaling behaviors are truly from the 
effects of interactions between the electrons and the gauge boson. 
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Systems 

Excitation 

Coulomb 

2D Dirac [HEU 

marginal q.p. 

iso., marginally irr. 

3D Dirac i .‘Ml 

marginal q.p. 

iso., marginally irr. 

3D Quadratic [8] 

no q.p. 

iso., relevant 

3D Anisotropic [10] 

q.p. 

aniso., irr. 

2D Anisotropic 

marginal q.p. 

aniso., marginal 


TABLE I. Comparison with the quantum criticalities in various semimetallic systems. Here the 
second column represents types of allowed excitation, “q.p.” is for quasi-particle. The third column 
represents characteristics of screened Coulomb interaction, “iso.” is for isotropic, “aniso.“ is for 
anisotropic, and “irr.” is for irrelevant. 

DISCUSSION AND CONCLUSION 


The presence of the novel fixed point implies that the electrons and gauge bosons are 
strongly correlated. At low energy, electrons and gauge bosons affect each other, so the 
Coulomb interaction mediated by the bosons becomes anisotropic and electrons receive back- 
reaction from the renormalized anisotropic Coulomb interaction. Thus, the Coulomb inter¬ 
action behaves differently from that of most critical systems where it enforces low-energy 
isotropy of electronic modesJH] HU EH El]- Also, notice that the ground state of our fixed 
point has marginally well-defined quasi-particles as those in graphene, which is in contrast 
to non-Fermi liquids with non-zero anomalous dimensions. In table [T] the comparison with 
other quantum criticality associated with topological phase transitions is summarized. 

Our novel quantum criticality can be experimentally tested in the systems such as VOo- 
Ti02 heterostructure. Near the critical point, optical conductivity shows anisotropy inher¬ 
ited from the electron band structure. Straightforward calculation with current operators 
(jx,j y ) = {^a x ^,^Ak y a y ^) gives 


a xx (Q) oc 




(T yy{Q) OC a/O, 


(13) 


upto logarithmic corrections from the Coulomb interaction (see the supplementary infor¬ 
mation VII). However, as shown in the previous section, the screening charge due to the 
charged impurity, which can be measured in principle by scanning tunneling microscopy 
(STM), shows qualitatively isotropic behaviors. Such discrepancy between the two exper- 
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iments is a smoking gun of the novel quantum criticality in addition to thermodynamics 
quantities such as specific heat. 

It is worth to mention that disorder scattering in the non-interacting electrons Q is 
relevant[35], so our results work better for cleaner samples. We expect that there will be an 
intriguing interplay between the anisotropic Coulomb interaction and impurity scattering at 
the novel critical point, which we leave for the future problem. 

In conclusion, we have investigated the quantum criticality of the anisotropic semimetal 
which can be thought as the critical point between topological insulators and Dirac semimetal 
in two spatial dimensions. At the low-energy limit, we found the novel fixed point out of the 
interplay between critical electron modes and the long-range - Coulomb interaction. The 
non-trivial anisotropic renormalization of the Coulomb interaction and the logarithmic cor¬ 
rections manifest at various physical quantities including screening charge when the impurity 
charge is introduced. Surprisingly we have shown that the scaling behavior of the screening 
charge in distance from the impurity is isotropic despite of the underlying anisotropic nature 
of the system. 

-Note added : After the completion of the paper, we became aware of the independent 
work by H. Isobe, B.-J. Yang, A. Chubukov, J. Schmalian, and N. Nagaosa [2SJ. Similarity 
and differences between our work and theirs are discussed in supplementary information. 

-Author Contributions : G.Y.C. performed most calculation and E.-G. M. came up with 
the original idea. Both authors analyzed data and wrote the manuscript. 
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SUPPLEMENTARY INFORMATION FOR “NOVEL QUANTUM CRITICAL¬ 
ITY IN TWO DIMENSIONAL TOPOLOGICAL PHASE TRANSITIONS” 

I. POLARIZATION BUBBLE 

In this supplemental method, we present the detailed calculation of the polarization, or 
the boson self-energy Il(fc) appearing in the main text. For convenience, we assume iV-copies 
of the electrons coupled to the gauge boson. If we are interested in the only one copy, i.e., 
physical limit, one can simply take N — 1 at the end of the calculation. 


A. Deivation of Polarization 


In the leading order in 1/N expansion, the boson self-energy is simply the two-leg one-loop 
diagram. The fermion propagator appearing in the loop diagram is 

1 


9f(k,w) = 


(14) 


-iu + H{k ) ’ 

where H{k ) is the free fermion Hamiltonian equation ( 6 ) in the main text. It is convenient 
to write the propagator in the following form 

1 ^ 1 


g f {k,uj) = 


—iu + H{k) 


£- 


s=± 


iu + E s {k ) 


P s (k), 


with the energy E s [k ) = sE{k) = s^/v 2 kl + A 2 ky and the projection P s {k ) 


Ps{k) — -(1 + s 


H(k K 

E{k))' 


Now the polarization is given by 


n (f»,Q) = ^ l Tr g f (k + P,u} + n)g f (k-P,u) 
Performing the contour integral over u j and the trace Tr [• ], we finally have 

= _ e2N t ^, P )( £ (* + S ) + g (*-i )) 

A ^+(E(k + l) + E(k~D) 

in which V(k,p) is the projection operator 

2A=1,2 + f - 2 ) 


2 ’ 


V(k,p) = 1 - 


E(k + *)E(k - |) 


(15) 


(16) 


(17) 


(18) 


( 19 ) 


where ei(fe) = vk x and 62 (k) = Ak 2 are the appropriate functions. 
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B. Calculation of Polarization 


To evaluate the function Il(fc) = Il(fc,u; = 0) in equation ([l8|), we use the change of 
variables as following: 


k x = p x x, and k y = 

which gives rise to the following expression. 


v\p x 


A 


-y, 


n (p) = 


e 2 N v\p. 


47t 2 u V A 

where f = \p y \yj^^ and f$(x,y) is the function 

1 


dx dy fs(x,y), 


( 20 ) 


( 21 ) 


Siip,y) = 


—oo J —oo 


* 2 -\+{y 2 -^) 2 


\J 0+)p 2 +(y+§) 4 \/ (a:—| ) 2 +(y—§ ) 4 . 


\j{x + \) 2 + (y + ^y + \J(x - \Y + {y - 

We introduce a function 

/ oo poo 

dx dy ft(x,y). 

oo J —OO 

Then it is easy to see F(£) = F(—£), i.e., F(£) = F(|f |). Now, we can see that 

F(|f| -> 0) -> c t , and F(|£| -> oo) -> Cy£, 


( 22 ) 


(23) 


(24) 


where c x ~ 2.7 and c y ~ 2.5. We plot F(£) in the Fig. pi From the numerical calculation, 
we can extract both the constants. With these in hand, we can write the polarization as 


n(p) = 


e 2 N F(|f|) e 2 N 


-\ P. 


Att 2 v' y |f| 
aN 

•VlPylG(lfl), 


47T 2 U 


|Py|G(|f|) 


where we have used G(|f |) = F ^' 1 and a = From equation ([24]), we can see that 


Id 

G(|f | —> 0) — y —, and G(|f | —> oo) —> c y , 


(25) 


(26) 


and G(|f|) > 0 for any f. 


II. SCREENING EFFECTS IN QUANTUM CRITICALITIES 

Here, we summarize bare and self-energy of Coulomb boson propagators in various quan¬ 
tum criticalities. The one-loop boson self-energy is obtained in the standard way. With 
different electron Hamiltonian for each case, we obtain the self-energy which is illustrated 
in Table HH 
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FIG. 3. Plot of F(£) 


Systems 

Bare 

Boson Self-energy 

2D Dirac [H CD] 

\q\ 

\q\ 

3D Dirac [3lj 

q 2 

q 2 log A 

3D Quadratic [S] 

Q 2 

\q\ 

3D Double Weyl[TQ j.T2] 

q 2 

q± log A + \q z \ 

3D Anisotropic [ID] 

q 2 

Agi 2 + \q z \% 

2D Anisotropic 

\q\ 

V\q~x T + Wy\ 


TABLE II. Bare and self-energy of boson propagators in various quantum criticalities. The second 
and thir column represent the bare and the one-loop self-energy of various quantum criticalities. 
A is the UV cutoff and only schematic behaviors are shown. Note that 3D quantum criticalities 
depend on the UV cutoff except the non-Fermi liquid quadratic band touching. 
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III. FERMION SELF-ENERGY AND RG EQUATIONS IN STRONG-COUPLING 
LIMIT 


Here we compute the fermion self-energy at the large-N limit to extract the RG equations. 
We start with the expression for the fermion self-energy with the corrected boson propagator 


£/(fc,w) 


r H(q + k) _ 1 _ 

2 J q E(q + k) \q\ + ^f-\q y \G(^ q ) 


(27) 


To extract out the correction to the bare parameters {v,A}, we expand the in k 

near k — 0 , 


-E f (k) = vk I a I I I + Ak 2 yi,. 


(28) 


Here {I x ,I y } are the integrals 


lx 87r 2 

f = A 

y 8tt 2 


i , A 2 q y l 

qx qy [v^ + ax?' 2 \q\ + <f\ qy m q y 

v 4 q A -5v 2 q 2 x A 2 q A 1 

Qx (ly [v 2 q 2 + AX]^ Iql + 'fbvlG&Y 


(29) 


Apparently, the integrals are divergent and needed to be regulated. We regulate the integrals 
by 


J dq x dq y -A 



dq. 


y> 


(30) 


where the UV and IR cutoffs {A, /i} are introduced. 
With the expression equation ([28]), we find 


Sv 


SA 


_ 1 T \ JE f (k,uj) 

Tr [cr°] L 5{k x ) 

vl x 


k ——>-0 


_ 1 t^\ y SE f (k,(jy 

Tr[<7°] L d(k 2 ) . 

Aly 


k —^0 . uj —^0 


(31) 
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A. Strong-coupling limit 


We first evaluate the integrals equation (|29|) when 


a —>■ oo. 


x 47 r 2 aN 


v 47 r 2 aN 


dq x 

dq x 


dq, 




”[« 2 9l + A 2 q$] 3 / 2 |<fc|G(£,) 

u 4 </J - 5u 2 </ 2 .4 2 </ 4 1 

dq v - - - 


( 32 ) 


Vql + AXr ' 1 te|G«„) ’ 

in which the cutoff scheme equation (|30[) is assumed. We demonstrate the detail of evaluations 
only for I x . 


where 


L = 


4ir 2 aN 

„2 


dq x / dq v 


A 2 qt 


[a 2 ? 2 + A 2 <ri] 2 / 2 |g„|G(£, 


dq x f dq y 


A 2 q: 


47 T 2 aNv J q x J \q y \ + ^]3/2 G(£ q ) 




e 3 


Qx J —oo 
oo 


= llog(-) / 
^ h Jo 
2J X , ,A 

= ^ log h’ 




(i+a 3/2 G«) 

e 3 i 

(i+a 3/2 c© 


J X 




e 3 


(l +f 4)3 /2GK )' 


(33) 


(34) 


The integral is well-defined in that it has no divergence in the integral because 


£ 3 


(i+e ) 3/2 g(o 
e 3 i 


ii i 

yy— oc — 
^ c x ^ 


(i + e ) 3/2 g to 


^e—r f 




(35) 


The integral can be evaluated numerically by using the numeric value of G(£), and we have 
obtained J x ~ 0.18. The numeric value is not important but it is important to remember 
J x > 0 because of G(|£|) > 0. 

Similarly, we can write I y as following. 


/ = A 

y N g V j ’ 


(36) 
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where J y is the definite integral 


J y — 


d( 1-5?* 


i o f(i+?r /2 G(o' 


(37) 


The integeral is well-defined because 

l-5£ 4 


(i + £ 4 ) 5/2 £G(0 

1 - 5^ 4 1 


-5 1 1 

~p ,— oc 77 

^°° s c y s 


-y 

1 1 

>o £ c x f £, 


1 

oc —. 
Or 


(i+e) 5/2 {G(o «-» ' ic x /i~ c*- (38) 

Hence there is no more IR divergence than log(A//i) appearing in equation ([36]) in the renor¬ 
malization of A in the strong-coupling limit. The numeric value of J y is obtained as 0.03. 

Hence, under the RG step, we find that the velocity v and the inverse mass A are renor¬ 
malized as follows. 


v v + Sv — v(l + log(A)), 

2 7 

A->A + 6A = A( l + ^log(A)), 

from which we deduce the RG flow equations, 

dv 2 J x 

~dl = ~N V ’ 
lLA 2J y 

dl N ’ 


(39) 


(40) 


which have been shown in the main text. 

Notice that the frequency dependence is trivial since only the instantaneous Coulomb 
propagator is used. If one incorporates frequency dependence, only the numeric numbers 
J x , J y are modified but not the structure of our calculation. The absence of the frequency 
dependence indicates no vertex correction by the Ward identity. 


B. Finite Coupling Constant 


Here we will discuss the fate of the finite coupling constant a which is now slightly away 
from the strong-coupling limit. In the prsence of the finite coupling constant a, the bare 
term ~ \q\ in the boson propagator 


9b(q,w) 


1 

g|~n(g) 


(41) 
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cannot be ignored. Hence the above calculations of a —>■ oo should be properly changed to 
see the scaling behaviors of the parameters {u, A}. 

Performing the straightforward calculation as the above, we find 


Sv 

v 


aK x (A,p), 


6A 

~A 


aK y (A,n ) 


(42) 


The integrals {K x ,K y } are 

_ /' A dqx , _ q+U 3 ' 2 _ 

1 J, q* J o V +^) 1/2 + fWo’ 

f A dq^ _ (i+^W 2 _ ^ 

w J, Qx Jo He+A{ qx )y/*+ { } 

where A(q x ) = ^ is the dimensionless number depending on the momentum q x . 

Before marching further into the details, we first make an assumption, A(q x ) < H(A) -C 1. 
The condition H(A) 1 should be interpreted as following. First of all, vA can be roughly 
thought as the bandwidth where the low-energy Hamiltonian equation (6) in the main text 
is the qualitatively correct description of the system. Note that A determines the curvature 
along k y and v determines the velocity along k x . Thus A(A) <g; 1 means that the velocity 
along k y is always much smaller than the velocity along k x within the energy E < vA, i.e., 
the dispersion along the k y direction is always much flatter than that along the k x direction. 
In other words, the spectrum is extremely anisotropic in that the equal energy contour is 
extended along k y and shrinked along k x . 

At this stage, the assumption A(A) <C 1 looks like an ad-hoc assumption to evaluate 
the integrals equation([43|). However, employeeing the momentum-shell RG approach at the 
weak-coupling limit a —> 0 (as present in the main text and the next supplemental method), 
we can show A 1 at the fixed point. 

With this condition in hand, we proceed to evaluate the integrals equation (|43| approxi¬ 
mately by taking A —y 0, 

K x = log ( — )J x (a) 

Ky = log (y)Jy(a) (44) 

/i 

where J x > 0 and J y are the constants depending on a whose detailed calculation can be 
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found below. From these, we find 


Hence we find that 


Sv 

v 

5A 

~A 


a log(—)J x (a), 
h 

alog(-)J 2/ (o;). 

/i 


jjogH 
log (A) 
^log(A) 
log (A) 


OtJxifA) 

aJy(a). 


( 45 ) 


(46) 


Thus the parameters {u, A} in the bare fermion Hamiltonian equation (6) in the main 
text receive the logarithmic corrections, and the coupling constant a decreases along the 
renormalization group process. 


1. Calculation of J x and J y 


Below we only present the detail for evaluating J x but it is straightforward to generalize 


to Jij 


J x = yyy j dq x dq. 
i 2 f A dq. 


A 2 ql 


8n 2 J ■■ ’■■■■>y \ v 2 q 2 + A 2 q^f/2 \ q \ + *K\q y \G(£ Q ) ’ 


•oo 

— I dq y 


A 2 q± 


4vr 2 u J y q x J_ ^ y 


a 

2 


r*A 


dq x dq v 


(1 + Sf ) 3/2 M + 

{“ 


- 4 1 


a 

M 


qx J- oo \q y I (i + £ 4 ) 3/2 i + | + °^G(f) 

= q / dq x r co di_e _i_ 

A q x Jo £ (1 + £ 4 ) 3/2 + dip _|_ e^G(0 

in which A(g x ) = ^q x and q y = We continue to evaluate this 


j f dq x 

J r = a - 


df e 


q x Jo 


= a 


£ ( X + ^ 3/2 y 1 + ^4 + °^G(0 

e 4 _1_ 

& Jo (1 + e) 3/2 y/e 2 + A(g x ) + !f£G(£) 


dq x , , c 
— / d£ 


(47) 


(48) 
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Now assuming A(q x ) < A(A) <C 1 which we will justify in the low-energy limit later, we can 
ignore the dependence on A in the integral to obtain 


J r ~ a 


dq x 




e 


q x J o s (1-K 4 ) 3/2 £ + ¥&(£) 


^r^Td, * 


q x J o s (l + a 3/2 i + ^G(0 


= a log(—) J x (q!) 
A 


where 


Jxipd) / d£ 


e 


1 0 S (l + e) 3/2 l + ^G(e q )' 

This integral is well-defined and hnite as 


e 3 


(1 + e 4 ) 3/2 1 + ^ G (^) 

e 3 i 


i 

£->°° £ 3 


(i+a 3/2 1 + 






Also J x (a) is always positive as G(£) > 0. 

On the other hand, we can evaluate J y similarly in A -C 1 and obtain 

A 


in which 


Jy{oi) 


J y = a log(—) J y (a), 
A 


<% 1-5? 1 


I o f (l+f‘) 5/2 l + !fG(f)' 
The integral is hnite and well-defined as 


l-5£ 4 


(l + £ 4 ) 5/2 l + 

1 - 5^ 4 1 


(l + e 4 ) 5/2 l + ^fG(0 


°C 77 

1 

oc —. 
£—■ > o c x 


(49) 


(50) 


(51) 


(52) 


(53) 


(54) 


Hence there is no more divergence in the renormalization in the inverse mass than oc log(-) 
in equation([52|. 
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IV. DETAILED CALCULATION OF MOMENTUM SHELL RG 


In this supplemental method, we will perform the one-loop RG by the momentum-shell 
method. We start with the following bare action S = present in the main 

text, 


S* = / ^ ku (-ioJ + vk x a x + Ak 2 a y )^ k)iJ , 


<S<P — 
Sip,<P - 


' k,uj 


' k,uj 


'VQl + ^ +«9*)|0(9,w)| 2 , 


' k,uj J q,Q 




(55) 


where rj is introduced to facilitate the engineering dimensional analysis of the anisotropic 
scaling along x and y in the bare boson propagator. The physical value of 77 , as mentioned 
in the main text, is 1 because the (bare) Coulomb interaction is isotropic. On the other 
hand, we have introduced k into the bare action to allow the possible anomalous screening 
term ~ \J\q x |, via the anomalous scaling dimension appearing in ~ q 2 term, under the RG 
steps. 

We start with the dimensional analysis. First of all, we set the engineering dimensions 
[x] = —zi, [y\ = —1 and [r] = — z 2 . Then it is straightforward to show 

[v}= z- z lf [A] = z- 2, [e 2 ] = z - 

W = 1 - Zi, [A] = 1, = (56) 


where A is the momentum cutoff in the momentum q y . Here the fields have the following 
dimensions. 


On the other hand, in terms of physical unit, we have 

M = M = A, [e 2 ] = 

H = !. [A] = p M = L 


The field operators have the following scaling dimensions. 


# = 


1 




i 

7z 


(57) 


(58) 


(59) 
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From the above dimensional analysis we have the following dimensionless constants, which 
will be used to parametrize the renormalization flow 

e 2 y/rj _ Arj A RA 2 A 3 y/fj 
2tt 2 v ’ v ’ ^ v 

Before proceeding to the details of the RG calculation, we first note that ij and e are pro¬ 
tected under the RG process due to the non-analytic structure of the propagator of the boson 
and the gauge invariance. Hence, the only non-trivial corrections are the renormalizations of 
{u, A} in the fermion propagator and k in the boson propagator. The renormalization can 
be deduced from the self-energies of the fermion and the boson, which we calculate below. 

We start with the fermion self-energy. 



—Ef(k,uj) = vk x a x h + Ak 2 a y I 2 , 


(61) 


such that 


h = ~ 
2 


h = ~r 


Aq 


ml + f + ml 

1 


- 3/2 ’ 

v 2 q 2 x + A 2 q* 
v 2 ql{v 2 q 2 x ~ 5 A 2 q A y ) 


ml + f + W x 


v 2 ql + A 2 q 4 


5/2 


where the integrals should be properly regularized as following. 


47T 2 ( 


dq y + 




dq T 


i —A 


r A 


dq x = — / dq y 


7T- 


dq x , 


(62) 


(63) 


as the integrands of equation (|62| are even under qt —» — q i: i = x,y. Here /j = Ae~ l ~ 
A(1 — l) with Z 1 is assumed by following the standard strategy of the momentum-shell 
calculations. 

With the assumption 4 < 1 which will be justified self-consistently in the end of the 
calculation, we can evalute the integrals straightforwardly as following. 


where 


h = al-Fibt 

h = alFiYj], 


UN 

fM - y ] 



dx 

dx 


1 1 

1 + yx 2 ( x 2 + l ) 3 / 2 ’ 
1 x 2 (x 2 — 5) 

1 + yx 2 (x 2 + l) 5 / 2 


(64) 


(65) 
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We remark that in 7 —>■ 0 limit (k —> 0), the inverse mass correction F 2 diverges logarithmi¬ 
cally in ~ logy. This indicates instability of the non-interacting fermion ground state with 
the bare Coulomb interaction without operator insertions. 

We now compute the polarization by the momentum-shell integral method 

-n(p) = plh, (66) 


in which 

e 2 f A 

i dky 

Now we have the renormalized action S' = S' p + S' 0 + S' p ^obtained by integrating out the 
high-energy modes as the following. We start with the fermion part S' p where 

d 3 x -0(<9 O + ie </>)-0 + J d 3 x 7p(v(l + aF 1 [y}l)'y 1 d x - A{1 + aF 2 [y}l)(d 2 ) 2 )'ip. ( 68 ) 

On the other hand, the boson part S\ is 

S* = f d ^ x l^\J 7 ^ + < ~ + ( 69 ) 

where 



dk 7 


v 2 A 2 ky 


a 


[v 2 k 2 + A 2 ky] b / 2 3 A 2 A 3 ^/rj/v 2 


(67) 


P=\~ 

37 


(70) 


On the other hand, the vertex correction vanishes and thus is the same as from 
equation([55|. 

From these, we can obtain the flow equations. 


d 

Jl 

d 

Jl 

d 

Jl 

d 

Jl 

d 

Jl 


log V — Z — Z\ + Jl 
log A = z - 2 + J 2 
log 77 = 1 — Z\ 
log K = ---h p 


log e 2 = z — 


l + ^i 


(71) 
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We now derive the following renormalization flow equations 

dot 2n r 1 

di = ““ 

^ | +7[-3 + 2a(F 2 [7] - Fifr])], 

^=/3(3- / 9)+a/3[F 1 [ 7 ]-2F 2 [ 7 ]]. 


(72) 


The fixed point of these equations is ( 0 , 7 ,/3) —> (0,0,3). We plot the numerical solutions 
also in Fig|4j This also implies that (ctFipy], aF 2 [y]) —> (0,0). With this fixed point, we 

>7 + $1 + 

d 


now derive the fixed point action S* = SL + S{ + 


^ + *>£*< 7 * + Ak 2 y a y )^k,u, 

J k,bJ 


s{ = 


' k,uj 


o 9 + 




S i,<t>= I I ie AH, 

J k,uj J q,£l 


(73) 


where we have used 7 = 1 for physical limit and (1 + /3/) ~ gj:e 3i —* y 7 1 | because 

i3=3 

e~ 21 is the scaling factor of q x . 

Now we justify 1 which has been used in the derivations of the fermion self-energy. 
We calculate 

— log A = — 1 + J 2 — Ji (74) 

and we can approximate J 2 = J\ ~ 0 near the fixed point. Then we clearly have 

— log A — 1 (75) 

which dictates that A —)■ 0 -C 1 at the fixed point. 


V. STABILITY OF NEW FIXED POINT 


In this supplemental method, we diagnose the stability and the scaling properties near 
the new intermediate fixed point. To determine the stability of the fixed point and the 
behavior near the fixed point more carefully, we now take equation([73| as the bare action 
and perform the momentum shell renormalization analysis again. Taking equation([73| as 
the bare action, we have 

M = * = l/z- 


(76) 



(a) a(l) 


(b) y(0 



(c) m 



FIG. 4. Plot of RG flows. Here l represents the RG time. We have taken the initial conditions 

q[0] = 7[0] = 0.1. 


Within the theory equation ([73]), we can introduce a dimensionless number 


hi = hi 



(77) 


which will determine the renormalization flow near the fixed point equation ([73]). Due to the 
non-analyticity and the gauge invariance, {e, 77 , n} are protected under the renormalization 
process. Then the only non-trivial renormalizations occur in {w, A}, coming from the fermion 
self-energy. We again evaluate the fermion self-energy 


= vk x a x I[ + Ak 2 y a y I' 2 


(78) 
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in which 


e 

I[ = - 
1 2 


6 

v — — 
lr > — ~ 


A 2 qt 


+ K\/\qx\ v 2 ql + A 2 cli 


3/2 


q l 

7j ' V IJ-I l/ 'ix 1 J A 'iy 

1 v 2 ql{v 2 q 2 x - 5A 2 q J) 


VQx + q ~w + *’ 2 ? 2 + 


5/2 


where the momentum integral J is regularized as in equation(|63|). 

After the straightforward algebra with A <C 1 near the fixed point, we find 

I[ = alfi[R] 

1 2 = alf 2 [R\ 


(79) 


(80) 


where 


fM = 

h\A = 


dx- 


Jq 1 + Ky/x (iC 2 + l) 3 / 2 
r°° i x 2 (x 2 — 5) 

J 0 dX l + Ry/^(x 2 +l) 5 / 2 ' 

From these, we can derive the flow equations 

^ logn = z-Z! + afi[R] 
Jl l °gA = z-2 + af 2 [R} 

jj log ?/ = 1- Z! 

d 1 

^ log ^ — 2 

d 9 l + ^i 

- loge =z- — 


(81) 


(82) 


which implies 

= -« 2 /i[«] 
n cv 

ji R= 2 ( 83 ) 

The attractive hxed point of these flow equations is (a, R) —> (0, R* ~ 0.35) such that 
f 2 [R*] = fi[R*]- The numerical solution can be found in Figj5j Hence we find that the 
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(a) a(l) 


(b) k(Z) 



(c) a m)\ 



FIG. 5. Plot of RG flows. Here l represents the RG time. We have taken the initial conditions 
q[0] = k[0] = 0.5. 


velocity and inverse mass at the energy scale /i are 




vi/j) & v(l + aCf log(—)) 

I 1 

A(fj) ~ A( 1 + aCf log(—)) 

/i 


(84) 


where Cf = A A*] ~ 0-8- This logarithmic divergence as fi —> 0 is reminiscent of graphene 


case. 


VI. DETAILED CALCULATION OF SCREENING CHARGE 

In this supplemental method, we show the detailed calculation of the screening charge 
dictated by equation (18) in the main text and equation (20) in the main text when a single 


31 








impurity charge Z is introduced at r — 0. We calculate the free and interacting cases at 
the level of the linear response theory. The details of the calculation here will closely follow 
Biswas, Son, and Sachdev[35]. 

2. Free Theory 

We first compute the screening charge of the non-interacting theory from equation (7) in 
the main text. The screening charge is given by 


pLM) = ^ B o(g) n o(<j) 


(85) 


in which D 0 (q ) = ^ and n 0 (g) is the one-loop polarization. The real space configuration 
of the screening charge, we perform the Fourier transformation 



( 86 ) 


We will be mainly interested in the directional behavior of the scaling in the integrated 
screening charges Q x (x) and Q y (y ) along x and y, which are obtained from 



— OO 


'OO 


d v pLd ( r ) 

dx p{ nd (r) 


Qv(y) = J 


'OO 


(87) 


—OO 


We now evaluate these screening charges. We start with 



( 88 ) 


Now from 



ZD o (q X} 0)U o (q x ,0) i 

1* = \i | 


z e 2 c x 1 
4tt 2 Vm4 y/\q^\ 9®=^’ 


(89) 
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we can take the derivative with respect to x on the both sides. After the straightforward 
calculation, we have 


Qx(.x) 


Zac T Pv 1 


8vr 2 V A yf\x\ 

Similarly, we have the following expression for the other direction 


Q v (y) = z (A)(o, 9j/)n 0 (o, q y ) 


qy ~\v\ 


By plugging 


Dq{0, 9j/)n 0 (0, q y ) = 


ec y 
4ir 2 v ’ 


we find that 


(90) 


(91) 


(92) 


Qy(y) X S(y). (93) 

Hence we see that the free theory has completely different scaling behaviors of the integrated 
screening charge along x and y directions, i.e., 

Q x (x) oc —^=, Q y (y) oc —Z5(y). (94) 

v x 

whose sign is the opposite of the impurity charge Z. 


3. Interacting Theory 


Now we use the non-trivial anisotropic fixed point equation (18) in the main text obtained 
from the interactions between the fermion and the boson. We will use the renormalized boson 
propagator dictated from equation (18) in the main text and the scaling behavior equation 
(20) in the main text to deduce the asymptotic behaviors of the integrated screening charges. 

PL(<?) ~ ZD vem (q)U 0 (q), 

D rem (q)~ 1 = \q\ + Ky/\q x |. (95) 


From this, we now calculate the integrated screening charge Q x (x). To calculate it, we 
need 


n 0 (^, o) 


e 2 c x 

4vr 2 v(q x ) 



(96) 
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in which {v(q), A(q)} are the running parameters at the scale q. From equation(j84|), we have 


n 0 (g*,o) = 


e 2 c a 




47 T 2 v(q x ) V A 0 
where {uo,H 0 } are the bare parameters. Hence we find 


Pind(O'x) Z 


e 2 c. 


v 0 


(97) 


(98) 


47T 2 Kv(q x )\ H 0 ’ 

where we have v(q x ) = Vo(l + aCf log(^T-)) with the 1 ~ E 0 , the bandwidth of the bare 
action equation (6) in the main text. From the relationship 


we find 




r / i 

/ Qx(x ")dx ~ Pind(0'x ] r)j 

1 0 fI 


Z(aC f ) 2 c 

2 /i[k*] V |x|(l + aC/log(|a;|/ro)) s 



oc 


M(1 + aCf log(|a;|/r 0 )) 2 ’ 
whose sign is the same as the impurity charge Z. 

From 


may) = 


C Cn, 


An 2 v(q y ) 


l9y|, 


we find 


Thus we see that 


Pind(0, q-y ) 


e 2 Zc y 
4:ir 2 v(q y )' 


Qy(y ) 


a 


Cy) 2 ^^ 


2/i[k*] |p|(l + oFj log(|7/|/r 0 )) 5 


oc 


IpK 1 + aC/ log(|7/|/r 0 )) 2 
whose sign is the same as the impurity charge Z. 


(99) 


( 100 ) 


( 101 ) 


( 102 ) 


(103) 


VII. DETAILED CALCULATION OF OPTICAL CONDUCTIVITY 


In this supplemental method, we compute the optical conductivity crjj(Q),j = x,y at the 
level of the linear response theory. The expression for the optical conductivity is given by 


= ^ Im n jj(° + i6 ) 


, 3 =x,y 


(104) 
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in which II^(il+i5) = —>■ fl-M<5) where the latter is obtained from the imaginary-time 

formalism. Hence we need to compute 


n^iJl) = / Tr Mjgf(k,uj + Q)Mjgf(k,Q) , 

J k,Ld 

np(a) — np{/3 ) 


E L^{ M i p e,c M iP^)z 7 


a,(3=±l 


iVL + aEt - (3 Et 


(105) 


in which Mj is the appropriate vertex along the j-direction and the factor P^ a = |^1 + 
a^J,a = ±1 comes from the Green’s function gf(k,u>). By performing summation over 
{a,/3}, it is straightforward to find 


+ 


iQ + 2E k —iVt + 2 E^ 


-1 


%(«) = x (- 

Here we perform the Wick rotation i£l —> H + id. 

n«w) = x ( n + u _ „ + u + 2E 


Now it is straightforward to take the imaginary part from this expression, 
sumption fl > 0, we have 


(106) 


-) (107) 

k / 

With the as- 


Im 


nj(n)] =>r lix{M j p- K+l M j p- K _ 1 )s(n-2E k \ 


(108) 


which can be evaluated by performing elementary integrals over k. By combining with 
equation( |104[ ), we finally obtain 

= 77 - 2 E„). (109) 

•JJz 

Now we evaluate this for M x = vea x and M y = Aek y a y to calculate cr xx (Q) and cr yy (Q) for 
this anisotropic semimetal. 

Plugging M x = vea x into equation( |109[ ), we have 

cr xx (H) = J dk x dk y Tx[a x P^ +l a x Pj: _^j x 5{E k - ^) (110) 


By the change of variables ( x = vk x and Q, = Ak y and then performing the trace, we have 


&XX (“ ^ ) 


e 2 v 


igkVaq 


d(*dCA\/Q + C* - f) X (l - 


( 111 ) 
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Performing the change of the variables ( x = and Q y = \/hh/, we hnd 


@ xx(S^ 1 ') 


e 2 v 


16n^/A \/fi 

with the constant I xx dehned as following 

V = f dxdy c5(\/l 2 + y i - ))(l + 


(112) 


4 2 

y — x 
y A + x 1 


2.47 


(113) 


We can plug M y = Aek y a y to calculate a yy (kV). After the straightforward calculation, we 
hnd 

e 2 \/A /— 

a yy (fl) = —^——IyyVQ, (114) 

where I yy is a constant 

4 y= I dxd V S ^ x2 + y4 ~lK 1 + y 4 + l 2 ) : 

«4.94 (115) 


In summary, we found anisotropic scalings in the optical conductivities along x and y 
directions 

cr xx (fl) oc l/VT), <T yy (fl) oc y/Q, (H6) 


VIII. INDEPENDENCE OF CUTOFF SCHEME AT QCP 


In this supplementary material, we show clearly that the RG how equations near the 
QCP in |V) i.e., Eq(8) in main text and Eq(69) of[Vj are cutoff scheme independent. To 


show this, we hrst note that the only non-trivial renormalization to the bare action (73) is 
from the fermion self-energy, and thus we only need to show the independence of the cutoff 
scheme in the renormalization of the velocity v and the mass A, which can be calculated 
from 


-D A Y>f(k,uj) = -A-^E f (k,u). (117) 

We introduce a soft cutoff function J(^) = ,/(—4) which is a smooth function (which is 
monotonically decreasing for q y > 0) satisfying 


J(0) = 1, J(oo) = 0. 


(118) 
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Then the fermion self-energy can be written as 

v /. n e 2 f H(q + k ) 

= — j 


qE{q + k) ^ q ~X^ 


(119) 


We will explicitly demonstrate that the renormalizations to the velocity v and mass A are 
independent of the detailed form of the function «/(■), i.e., the flow equations near the QCP 
(see |V|) , 

~jj log v = z-Z! + afi[R\ 

~jj log A = z-2 + af 2 [R} 


( 120 ) 


with 


/i[«] = 
/a[«] = 


dx- 


l 0 1 + Ky/x ( X 2 + l) 3 / 2 

, 1 x 2 (x 2 — 5) 

ax- 


( 121 ) 


/ 0 1 + Ryfx ( X 2 + l) 5 / 2 ’ 

are independent of the precise form of the cutoff function J(-). Note that the flow equations 
in 0 are originally obtained by using the hard cutoff function. 

A. Correction to the velocity v: we first calculate the correction to the velocity. We 
only need to calculate 




H(q + k) 


9b(q)J 2 (^-)i 


q E (q + k)* UK ' 1 ''' v A 


with k = (k x , 0). Expanding in to 0(k x ), we find 

-E f (k,uj) = vk x a x x I x , 


where 


Ix= r 2 / dQx dqy 


A 2 qt 


J 2 (X 


[v 2 ql + A 2 q*\ 3/ 2 q y + K^Jq x A 


We now perform the derivative D\ on the self-energy to find 


-Z>aP/(*,w)] = vk x o x x D A [4]. 


( 122 ) 


(123) 


(124) 


(125) 


Hence we calculate 


e 2 f 

Da[I x ] = ^ J d 9x J dq y [ 


A 2 qt 


[v 2 q 2 x + A 2 g4]3/2 Qy + 


Da[J 2 ( f)]). (126) 
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To evaluate the integral, we first perform the change of the variable y = % to find 


e 2 f°° ( 

D A [I X \ = -—Ax j dq x J dy[ 


A 2 AV 


[v 2 q 2 + A 2 A 4 ?/ 4 ] 3 / 2 Ay + K^Jqx dy 


vAAv)\ 


Next we scale out q x —* q x y to find 

.2 roc , A 2 A 4 


D a [I x ] = 


2 vr 2 


A 


dq x 


i o 


[v 2 q 2 x + A 2 A 4 ] 3 / 2 A + J 0 dy dy ^ 


where the integral over y can be performed analyticall to find 

.2 roc , A 2 A 4 i 


e 2 

£>A I / J = l d <h V [,,2,2 + 212A4JS/2 A + K ,J, h p 


(127) 


(128) 


(129) 


which is nothing but the momentum-shell integral for the renormalization of the velocity v 
obtained in ©of0 he., I[ in Eq.([79|). Notice that the expression is manifestly independent 
of the precise form of «/(•). After successive change of variables as in[Vj we can show that 


the same integral fi[n] in Eq.(120) appears in the renormalization of the velocity v. 

B. Correction to the mass A: we next calculate the correction to the mass A. After 


the straightforward series expansion of (122) for k = (0 ,k y ) to 0(k 2 ) and then taking D Al 
we find 


-r> A p/(fc, w)] = Ak 2 y(K , + k 2 + k 3 ), 


(130) 


where {K\, K 2} K 3 } are the following integrals to be evaluated. 
- 2 /• ' 1 


K, = (a 7(|)7"(|) + |j'(|)j"(|) + |j(|)J'"(|)) ), (131) 


in which f(q) = 


Ml 


\Jv 2 q£+A 2 q* 

f"(x),f'"(x). For /l 2 , we have 


and we have used the notation f'(x) = 4r.f{x ) and similarly for 


*2 = £ ( A Mf )/(f) + ,„A[/(f)]’ + V"(f)) ), (132) 


and, for Ii 3 , we have 


K e2 f 9b{q ) a T( q y . V( q y . 

K ^2A.l q — f 


(133) 


Performing the change of the variable q y = Ay and then scaling q x —y q x y , we can show that 

d 


K\ oc J dyA-[y 2 J(y)J"(y)}=y 2 J(y)J"(y) 
r°° d 

K 2 (x J dy—[yJ(y)J'(y)\ = y 2 J{y)J\y) 


= 0, 


= 0 . 


(134) 


38 















On the other hand, we find 


K, = 


2vr 2 


A 


dq x 




5A 2 A 4 ] 


x 


'0 


'0 


+ A [v 2 q 2 + A 2 A 4 ] 5 / 2 
By performing the integral over y explicitly, we find 

- 2 '■°° 1 v 2 q 2 [v 2 q 2 — 5A 2 A 4 ] 


d 4y W 


(135) 


K 3 = -^A 

3 2t r 2 


dq x 


(136) 


' o 


+ A [v 2 q 2 + A 2 A 4 ] 5 / 2 

which is nothing but the momentum-shell integral for the renormalization of the mass A 
obtained in ®of0 he.,/' in Eq.([79|). Notice that the expression is manifestly independent 
of the precise form of «/(•). After successive change of variables as in[Vj we can show that 


the same integral / 2 [k] in Eq.(120) appears in the renormalization of the mass A. 


IX. WARD IDENTITY AND NON-RENORMALIZATION OF ELECTRIC CHARGE 


Here, we discuss the Ward identity and its consequences in two dimensional systems with 
the Coulomb interaction, and we also discuss difference between our work and Isobe et. 
al. [3E] • The effective action after redefinition e<j) —> 0 is 

S=( + i(j)) + 'H(-iV))il> f j^\4>(q,uj)\ 2 + S ins , 

J £C,T J q,UJ 

where S ms represents an additional inserted operator. For the discussion of the Ward iden¬ 
tity, the inserted operator part is not important but it becomes crucial to investigate the 
stability of the non-interacting fixed point. Note that the action has gauge-invariance under 

(j) —> (j) + d r a(T), -0 — > ' 0 e_ * Q( '" r ^ 


After incorporating quantum fluctuation, the effective action becomes 

S —> j{\ + + (1 + S v ) 

Notice that non-analytic dependence of the A\q\ term in q prohibits its renormalization 
under integrating out the high-energy modes. Here we did not explicitly write out the 
renormalized Hamiltonian H(—iV) as it is not important in the discussion of the Ward 
identity. 

It is clear that the gauge invariance gives the Ward identity, 


itfi'ijj'xl) + 


+ f H|0|2 + s„ 


= &v 


39 







Moreover the non-renormalizable e 2 gives that the correction terms are absorbed by redefin¬ 
ing -0, 0 = A y / Z^0 r with Z0 1 = 1 + S u , and parameters in the Hamiltonian H(-i'V) such as 
velocity and effective mass. 

Therefore, in our calculation, the electric charge is always non-renormalizable but eZ 0 
receives correction even though it appears in higher loop calculation with the bare Coulomb 
interaction. 

This is one main difference between work by Isobe et. al. |36j and ours. In their work, 
they claim electric charge (g in their notation) is renormalized but gZ ( eZ^ in our notation) 
is not. 

Another related difference between our work here and theirs[36] is their emphasis on the 
wavefunction renormalization Z — (1 + However, in the strong-coupling limit 

aN —> oo, it is not difficult to see that contains the log-squared divergence (see below 

for the detail), which originates from the poor screening of the Coulomb interaction at finite 
frequency with zero momentum. In this limit, the only effect of the Coulomb interaction 
is to rotate the phase of the fermions as well explained in Son [IT] and thus should not 
be taken as physical. Such log-squared dependence is also present in Isobe et.al. [36] in 
the limit aN —> oo and taken to calculate the fermion wavefunction renormalization Z. 
They interpret the wavefunction renormalization as single-particle residue. But it cannot 
be adiabatically connected to one of the strong coupling fixed point due to the additional 
divergence. 

One can see another difference from insertion operators S ms . As shown above, the pertur¬ 
bative RG with the bare Coulomb potential without the inserted operator gives log-squared 
divergence in the inverse effective mass correction. This log-squared correction indicates 
necessity of additional operators. In our calculation, we insert the operator 

S, n , = \ /(V*)’, 

which is natural in the Wilsonian RG, while they[32] insert an approximated polarization 
function even for weak coupling analysis, 

SL= f ^ n (g,a;)|0(g,a;)| 2 . 

We believe that the insertion of an approximated polarization function is questionable in the 
weak coupling limit since it is not controlled by the d- factor, which disappeared in their 
weak coupling analysis. 
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Note that the insertion of the approximated polarization which is contributed by the 
fermion modes from all energy scales is valid in strong coupling limit since it is well controlled 
by the factor. In the current problem, the strong coupling limit is indeed controlled, and 
all physical quantities receives corrections which are well captured in fsobe et. ah and 
ours with qualitative difference from different approximation. However, as shown in the main 
text, the strong coupling fixed point becomes immediately unstable due to the anomalous 
dimension of the velocity, thus the validity of the large -Nf expansion is not guaranteed in 
the infrared limit. 

We also believe that unusual cutoff (scale) dependence of physical quantities in their 
analysis (for example, log a in RG equations) indicates their RG scheme is not conventional 
and it cannot be perturbatively obtained from the non-interacting theory. On contrary, in 
our analysis, all physical quantities receive logarithmic cutoff (scale) dependence as ones of 
quantum criticality described by renormalizable held theories. 

Apart from the differences, we find the ground state and its excitation in the IR limit are 
marginally stable, so-called marginal Fermi liquid, which is also obtained by Isobe et. al. 
[36] in spite of different renormalization of the inverse mass. 


A. Log-squared Divergence in Fermion Self-Energy 

Here we present the log-squared divergence in the strong-coupling limit, aNf —> oo, where 
g^ 1 (q,Q) —>■ — n(qr, 0) by keeping the dependence in f2. 

Such log-squared divergence has been observed in the graphene case as discussed by 
Son[[H]. The fermion self-energy in Son is given by 

£(p,Po) = s o7oPo + £17 ■ P, 

where 

_ 8 _ f d?qdq 0 ql-q 2 

0 N f J (2tt) 3 (g 2 )3/2|q|2 

with q = | (g 0 , q) |, the size of the three-component vector (g 0 , q), which is not to be confused 
with the size of the spatial two-component vector \q\ = \(q x , q y )\. Naively, this expression 
is single-log divergent, but it is in fact log-squared divergent. The log-squared divergence is 


(137) 


(138) 
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apparent in the limit ^ —>■ 0 because Eq.(138) can be approximated 


£n OC 


d 2 qdq 0 


<7o 2 


(2vr)3 (g 0 2 )3/2 \ q \2 


oc 


dqo 

Qo 


x 


d 2 q 


(139) 


which has the two sources of the logarithmic divergence, one from the integral over go and 
the other from the integral over q. 

We expect the similar log-squared divergence to be present if we keep the dependence 
of the polarization II(g, hi) in hi. To exhibit this explicitly, we reconsider the fermion self¬ 
energy and the boson polarization in the strong-coupling limit a —> oc. Here we show the 
log-sqaured divergence in E 0 where 


£o = --Tr[E / (fc,u/)]. 


(140) 


Here we work in the unit v = A = 1 for convenience and the clarity of the discussion. We 
start from the expression 


£n OC 


n 2 ~E 2 q 




(141) 


where E q = (g 2 + gH 1 / 2 . In the limit Q E q , we have 


£n OC 


,n ™ n( 9 ,si) K L n x J„ + tt ,,j v /p- 


(142) 


where a x ~ 0.66 and a v ~ 0.75 and we have used the following form of the polarization 


n(q,SJ)cc^ a ,-^= + a,4Vm 


(143) 


which is asymptotically correct as far as O > i? 9 . Now in Eq.(142), we can perform the 
change of the variables x = ^j 4 , q x and y = ^Oy\il\ x Aq y to find 


r dn 

E 0 oc / - x 


dxdy 
x 2 + y 2 ’ 


(144) 


which clearly exhibits the log-sqaured divergence. 

Now we show how we obtained the asymptotic expression of the boson self-energy 


Eq.(143). We start with the expression for the boson self-energy Eq.(18). By performing 
the following change of the variables k x = |f2|x, k y = |Q| l ^ 2 y, X — |% and Y = 


1y 

a yw\ 
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(remember that we took v — A — 1 for this subsection) and expanding hi to the lowest 
orders in X and Y, we find 


n (q, n) « \/\n\\a x X 2 + a y Y 2 


1 


: V^T 


+ ^VWlQy ■ 


(145) 
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